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Abstract 

We describe a nonperturbative calculation of the spectrum of SU(2) 
Yang-Mills theory based on a Haniiltonian formulation. Our approach 
exploits gauge invariant variables similar to those used in nuclear 
physics to describe collective motion in nuclei. 
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Our understanding of the low energy behaviour of QCD depends crucially 
on the developpment of nonperturbative methods. The variational gaussian 
approximation which has been successfuU in quantum mechanics and in scalar 
field theory remains difficult to apply for a nonabelian gauge theory because 
of the requirement to satisfy the Gauss law constraint. The necessity to 
maintain gauge invariance is however an essential ingredient in practical cal- 
culations. In order to study the configurations which contribute significantly 
to the nonperturbative ground state, we have to take into account the gauge 
invariant functional measure. The nontrivial gauge invariant volume element 
will induce a centrifugal effect and the boundary conditions on the wave func- 
tional will differ drastically from those in perturbative calculations. This has 
important phenomelogical consequences such as the occurence of nonvanish- 
ing condensates describing the vacuum state and the existence of a mass gap 

i- 

In this paper, we use a representation of the SU(2) vector potential which 
separates explicitly the gauge degrees of freedom 0, ^. In this representa- 
tion, the Gauss law appears as a local constraint. The hamiltonian however 
becomes non-local. A derivative expansion (or strong coupling expansion) 
valid in the nonperturbative domain and for slowly-varying fields allows one 
to write expliciltly the first few terms of an effective hamiltonian. One can 
thus obtain approximate solutions for the dynamics which are gauge invari- 
ant. 

In section 2, we introduce more appropriate gauge invariant variables 
p(x),/?(x),7(x). These variables describe the field configuration in an in- 
trinsic frame and they can be interpreted as "density" and "deformation" 
variables. They are analogous to the collective variables introduced by Bohr 
and Mottelson to describe the dynamics of deformed nuclei 0, Q, the pure 
gauge degree of freedom corresponding to intrinsic nucleon coordinates. The 
collective gauge invariant variables p, (3, 7 allow one to describe deformed so- 



lutions of the dynamical equation without breaking the local gauge symmetry 
of the hamiltonian. 

In section 3 and 4, we work in the strong coupling limit which corre- 
sponds to the constant field approximation. We give the expression of the 
hamiltonian in terms of the gauge invariant variables p, (3, 7. We investigate 
the properties of the ground state in the SU(2) gauge theory and we obtain 
predictions for the lowest masses of the color singlet bound states of gluons 
called glueballs. The minimum of the energy corresponds to a configuration 
with axial symmetry (7 = 0) and a strong deformation (/3 near 1). In this 
formalism, the lowest glueball states are interpreted as rotational levels and 
vibrational collective levels. In the last section, we discuss our results for the 
behaviour of the wavefunctional. 

1 The polar representation 

In the hamiltonian formalism, we choose the temporal gauge A^ = 0. The 
hamiltonian reads : 

H = - j d?xtr (e^ + B^) . (1) 

For the SU(2) color group, the vector potential Aia (where i is a space index 
and a is a color index ) is a 3 x 3 matrix. The polar representation is given 
by[ii: 

where A„,n = 1,2,3 are three numbers, f{Oi) and /i(0a) are orthogonal ma- 
trices parametrized by two sets of three Euler angles 6i and (j)a and g is the 
bare coupling constant. In terms of the 3x3 spin 1 matrices {Si)jk = —ieijk, 

f = exp {-iOiSz) exp {-i92Sy) exp {-1638^) . (3) 

The matrix / describes a rotation in ordinary space and the matrix h a 
rotation in color space. In general, A„, 6i and (pa are space-dependent. Under 



a local gauge transformation, the matrix h is simply rotated while / and 
A„ remain unchanged. Therefore, among the nine variables Aia-, six gauge 
invariant variables A„ and Oi are explicitly separated from the three gauge 
degrees of freedom </)„. 

In these new variables, the Gauss law operator 

G\x) = - D^E,^ , (4) 

becomes a local operator : 

G\x) = J\x) , (5) 

where J^ are the cartesian components of the color angular momentum in 
the laboratory frame expressed in terms of the Euler angles 0^ : 

Ji = -«<-cos0iCot02^— -sm0i— — + ^— — — — ^ , (6) 

( acpi d<p2 sm 02 <^03 J 

with similar formulae for J2 and J3 [^. For a color singlet state |\1/ > 

G\x)\^ >= J\x)\^ >= Q . (7) 

Therefore, the wave functional for a color singlet state depends only on the 
six gauge invariant variables : \l/(A„,(x), ^j(x)). The functional integration 
measure becomes 

\[DAt{x) = n l^lix) - Xl{x)\ YlDKi^) dfdh , (8) 

i,a n>m p 

where df = sin ^2 d6id62d6^, and dh = sin 02 (i0i(i02(i03. 

In terms of the new gauge invariant variables, the hamiltonian becomes 
non local. The expression for the potential energy ^i?^ is given by : 

V = ^{XlXl + XlXl + XlXl) + W (c^rl w^f , (9) 

where Wn are the three orthogonal vectors : w)^ = finXn- 



The kinetic energy is a non local function of the gauge invariant quantities 
and their derivatives. However, we can use a derivative (or strong coupling) 
expansion to write explicitly the first few terms. This expansion will be valid 
in the nonperturbative domain and for slowly-varying fields. In the following, 
we will consider only the lowest order in 1/ g"^, which corresponds also to 
the constant field approximation. At this order, we keep only the constant 
term in (Q) , which is the non-abelian term. This approach is therefore not 
adequate to describe the abelian limit. 

For color singlet states, the kinetic energy is given by [S] 



^- -^m\Yl K\ x\ +2 E T^ Ti"h"7^~^ 



2-^ i„=l (^An^An n>m A„ — A^ \ 5A„, SXm^ 

1 \2 _|_ \2 

4:L\^ (A2 -A2)2 "-'= ^ ' 

k,n,m \ n mJ 
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(10) 



We have introduced a length scale L, the total volume being L'^. The opera- 
tors Lfc are the components of the angular momentum in the proper frame. 



As an example |jTO[ : 

. f COS03 d . r, d n /i"9l ,^^, 

In lowest order in l/f?^, the angular momentum in the laboratory Mi = fikL^ 
is equal to the spin density : M(x) = A"-{x) x E°-{x). 

It can be also usefull to write the expression of the gauge invariant oper- 
ator B"-.E°'. For color singlet states, it is given by : 



-D .H/ i g ^npq ^mpq ^p^ 



"'"ax ^ T2 T^-^m-t^fc 



;i2) 



2 The gauge invariant collective coordinates 

p{x),l3{x),-i{x) 

States of zero angular momentum depend only on three gauge invariant and 
rotational invariant quantities which can be choosen as [|I^ : 

B^ = Btix)Btix) , (13) 

B.y = Btix)ytix) , (14) 

and 

y' = yt{x)yt{x) , (15) 

where yi{x) = e°''"^eikiBl{x)Bf{x). For constant fields, they are related to the 
gauge invariant variables A„ according to the formulae : 



5^ = / [XiXi + XiXl + A?Ay , (16) 

y' = Ag^XlXl (A? + A^ + A^) , (17) 

B.y = Gg'XlXlXl . (18) 

These expresions, as well as the functional integration measure (^, are in- 
variant under a permutation of the A„ and a simultaneous change of sign of 
two A„. 

In the following, we will use three gauge invariant variables p, j3, 7 defined 
as : 

p' = Xl + Xl + Xl, (19) 

Xl-Xl = ^p'(3sm'y , (20) 

2Xl~{xl + Xl)=2p^(3cos^ , (21) 

i.e., 

A?-^ = ^//5cos(7 + y) . (22) 



A^-^ = ^P^/5cos(7-|) . (23) 

Al-^ = ^p2/5cos7. (24) 

where < /? < 1 and, from symmetry properties, we can restrict the angle 7 
between and vr/S. The variable p has the dimension of 1/L. 

The vibrational part of the wave function ( i. e. the part of the wave 
function which is independent of the three Euler angles 9i) is a function of 
the three gauge invariant collective coordinates p, /5, 7 : \1/ (p, /5, 7). The vari- 
ables p, /?, 7 are interpreted as the analogous of the density and deformation 
variables used in the collective model of Bohr in Nuclear Physics [Q. They 
remind also the derivation of the collective hamiltonian for a system of N 
particules by introducing 3 N — 9 Euler angles and six collective coordinates 
which describe the shape, the dimension and the orientation of the system 
[§]. Since we start from nine degrees of freedom A^ for the Yang-Mills system, 
the description in terms of the gauge invariant variables corresponds to the 
treatment of the four-body problem ^j. The variables p, /3,7 will describe 
the gluon configuration in the intrinsic frame. The p vibrations correspond to 
density vibrations of monopole character. The j3 and 7 vibrations correspond 
to quadrupole oscillations. In general, there will be a coupling between the 
oscillations of the density p and those of the deformations /5, 7. We will show 
that p, (3 and 7 are convenient coordinates to perform practical calculations. 
Furthermore, they give a physical insight in the structure of the vacuum state 
and the lowest excited states. 

In the /3, 7 plane, the axis (3 = corresponds to a "spherical" field config- 
uration : Ai = A2 = A3. The axis 7 = corresponds to an "axial symmetric" 
field configuration : 

\l = \l=P^{l-P') , (25) 
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A^ = ^ (1 + 2(5) . (26) 



The point /5 = 1,7 = corresponds to two vanishing A's : Ai = A2 = 0, and 
A3 = p^. We point out that this configuration does not describe an abehan 
type field. Indeed, in the strong couphng approximation which corresponds to 
the constant potential approximation, we are not in a position to investigate 
the abelian limit. 

An arbitrary point in the /3,7 plane corresponds to a "triaxial" field 
configuration. 

The gauge invariant quantities in eq. ([T6| - |l8D become functions of p^,/3^ 
and /5^ cos 87 : 

|^ = /y(l-/5^), (27) 

i?-Z/ = /^p'Q-/5' + ^/5'cos37) , (28) 

and the expression for the volume element of the vibrational coordinates 

dT=\X\- A^l |A2 - A^l |A2 - A^l dAi d\2 rfAg , (29) 

becomes : 

dr = —^ p^p'^l sin 87! | det M\ dp dp d-f , (30) 

8v 8 

where the jacobian is given by 

^^^^^^D^Wa)^ 1 

D{pP^) Gv^AiAsAg^"^ ^ ^ 

A singular point occurs when one of the A„ vanishes or when two A's are 
equal. This will imply suitable boundary conditions on the gauge invariant 
wave function \E' (p, /?, 7). Let us now use the following rescaling : 

$ (p, (3, 7) = p''^' TT-VWJ^ ^ ^P^ ^' ^) • (^2) 

IA1A2A3I 

The wave function $ is normalized according to the integration measure 

/3^ I sin 87 1 dp dp dj . (33) 



3 Strong coupling expansion 

In terms of the new gauge invariant variables p, /?, 7 and after the rescahng 
(|32[), the vibrational hamiltonian (i.e. the terms in equations (^) and (p!0|) 
independent of the angular momentum operators Lk ) is equal to T^b + V 
with 

9p2 



Tvib - 2L3 \ p. 2 



+ -| (1 + /3cos37- 2/52) ^ + -^ (1-/3 cos37)^ 

+ 2^- o (3 COS 37 - 3/3 + 2/3 sin^ 37)!- 1 , 
p^/?^ sm 37 ^ ^ t^T I 



(34) 



and 



^-(p,A7)^L3|p^(l-/3^) + ^ + ^ ^^^-^^3^^^3^ .(35) 



^' 4 /. .2\ , 99 , 3 1-/3 ^ 

i-/32 + |/ 

We wish to stress that the last two terms in the expression of the poten- 
tial energy (^) and in particular the 7-dependence arise from the nontrivial 
gauge invariant factor in the integration measure (|30|) . This potential di- 
verges for /? = 1,7 = 0. Therefore, the wave function has to vanish at this 
point. For 7 = and /3 = 1 — e, we have, in the limit of small e : 

vl/(p,/3,7) = v^V'^(/''/^'^) • (36) 

In the limit of small deformations /3 <^ 1, the vibrational kinetic energy 
can be written as : 

r Q2 



rp ]_ 

-'■ vib 2L^ 



dp^ 



p2 



]3^di3[^d(3y /32sin37 d^ V "" ^d^^ 



(37) 
We thus obtain in this limit the expression of Bohr and Mottelson for the 
vibrational kinetic energy. 
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The next step is to find the equihbrium shape by minimizing the po- 
tential energy V (p, P, 7) with respect to p, P, 7. For this purpose, we have 
calculated the derivatives of the potential and the general expression for the 
stability matrix. 

1 /'i 

Spherical configurations {(3 = and L'^pl = g~'^^^ [y) ) correspond to 
maxima with an energy equal to : V' (po, /^o, 7 = 0) = 5.9 g'^^^/L. 

We have found a minimum at 7 = 0, /3 = /5, p = p, where (3 and p are 
solutions of the following equations : 



The two functions given by the right hand sides of eq. (^) and of eq. (p^) 
intersect at (3 very near 1. By writing (3 = 1 — e and by keeping the lowest 
order in e, we find at the minimum 

and 

T2-2 /3\2/3 1 /99\2/3 1 

Therefore the minimum of V corresponds to a strongly deformed field con- 
figuration with axial symmetry. We obtain : 

V' (p-^ /3 = 1 - 6, 7 = 0) = //^ ^ (^'^') /L . (42) 

The corresponding value of the magnetic field is : 

B' = l9' -e-p' (43) 

or 

L^B^ = ^^y\m = gL-p . (44) 
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This minimum associated to a nonvanishing value of the magnetic field cannot 
be described by a perturbative approach which constructs a state around 

In lowest order in e or in 1/e, the eigenvalues of the stability matrix at 
the minimum are 

S= = 3© 9"'^, (47) 

where we have used dimensionless quantities. To this order in e, the eigenvec- 
tors of the stability matrix are given by the directions of 6p, 6e and 6'~f = rj. 
Therefore, the quadratic expansion of the potential energy arround the min- 
imum yields : 

V\p, (3, 7) = V'ip, e, 0) + \LW^i6pr + hcu'^idef + hco^ . (48) 

4 Solutions of the rotation-vibration hamil- 
tonian 

In analogy with the procedure used in nuclear physics, it will be more con- 
venient to redefine the wave function according to : 

I sin S'vl-^/^ 

^(p,A7) = |sin37r/^$(p,A7)=p''/'^-Vn72*(P'A7) • (49) 

IA1A2A3I 

After this rescaling, the differential operator T^h is transformed into T^b + 
Vadd- Near the minimum, p = p + 6p, P = 1 ~ e — 6e, 'j = 1], we have : 
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and 



^^'' -~u?J^^- ^^^^ 



In T„jb, we have neglected terms proportional to r] -^. 



The rotational kinetic energy in eq. ([10| ) can be written as : 



T, 



rot 



t2 t2 7-2 

J\ J2 Jz 



(52) 



2L3 

The rotational and the vibrational parts of the energy are coupled to each 
other due to the p, /5, 7 dependence of the moments of inertia : 

.2 sin^(7 + 27r/3) 
Ji = 2,p IJ — , . , (53) 

1 — pcos(7 + 27r/3) 

^ 9^2/,2 sin^ (7 + 7r/3) 

^^ = '^^ l + /5cos(7 + vr/3)' (^^^ 

:^3 = 2p^/5^ ^ '^"'^ . (55) 

1 — p cos 7 

The collective hamiltonian in terms of the Euler angles and the variables 
p, /3, 7 is given by 

H = Tr,t + %,^ + VaM + V'{p,(5,^) . (56) 

The wave functions are normalized according to the volume element 

j3^ dp d(3 d'-f dVt where dVt = sin62 dOi d92 dO^ is the volume element of the 

Euler angles. 

In order to obtain approximate solutions for the dynamics, we expand 



the potential energy arround the axial minimum according to eq. (48). Near 
the minimum, the moments of inertia are given by : 

Ji = ^^2 ^ p' , (57) 

and 

2 



J3 = ^ r/^ (58) 
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This yields for the kinetic rotational energy : 

The last term in T^ot is a rotation-vibration interaction. The presence of 
this centrifugal barrier has a deep connection with the gauge invariance. It 
arises because when two A's are equal we can not define the rotation angles 
in the polar representation for Aia- As a consequence, there is in this case 
no dynamics associated to these degrees of freedom. 

By keeping Ji and ^2 constants, we will neglect further rotation- vibration 
interaction terms. For small angular momentum, we can treat them in per- 
turbation. We have checked that they are indeed small, which is expected 
since the deformation is large. 

The eigenfunctions of H are of the form |T^ : 

1 /9 

if {6p, 6e, 7], Oi) = ( 16^2^1+3^0) ) Rn, (Sp) Qnp (5e) Xn, (^) 

X {d'mk m + (-1)' d'm-k m) , 

(60) 
where Rn^ are the solutions of the harmonic oscillator of frequence Up. For 
the functions Qnp {5e) and Xn iv)y we obtain the following expressions : 

gn, {5e) = A exp (^-^ (Se + e)') /7„^ (-^) , (61) 

and 

Xliv) = B \v\y' r^^l^ exp [-\r^^ ,F, (^-n„ ^ + 1, Ar^^) , (62) 

where Hn are the Hermite polynomials, iFi is the hypergeometric function 
and A = 3/2e^. A and B are normalization constants. 

We thus see from eq. (^) that dynamical effects arising from the collec- 
tive kinetic energy shift the center of the gaussian away from the minimum 
of the potential energy at /? = 1 — e to the point P = 1. This point was 



already known to be the minimum of the magnetic energy eq. (WA) - However 
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an additional information we have gained from eq. (^) is that the boundary 
condition on the wave function is precisely that it must vanish at /3 = 1. This 
implies that only odd values of n^ are acceptable. A second boundary condi- 
tion on the wave function is that it vanishes at 1]"^ = e^/3. For jS near 1 and 
7 small, we have from eq. ([T8|) and eq. (^Sf) : (Ai A2 A3) = p^ (e^ — Sr^^) /9. 
From the expression for Xn iv) 5 we see that the wave function is indeed neg- 
ligible at the boundary r]"^ ^ e^/3. 

The functions (|60|) describe states of positive parity, K being restricted 
to positive even integers and / = 0, 1, 2 for K = and I = K,K + 1... for 
K ^0. 

The eigenvalues of the energy are given by : 

EiKupupn^ = (up + - ) -ujp +[np + -) Ep+ hn^ + -\K\ + 1 ) E^ 



2L^p 



^(/(/ + l)-ir^)i + I/' 



2 L ' ■ V ^ 47 ^ V '2' 
1 

f 



(63) 



0^0x2/3 1 



where 



and cUp is given by 

The p vibration is therefore softer than the (3 and 7 vibrations. Note that 
the zero point vibrational energy for the (3 vibration is 7/4 Ep instead of the 
usual value 1/2. 

In the strong coupling approximation, the energy is proportional to g'^^^. 
By using e = 2/33, one finds : 

u;p=l.2Ag^/\ (65) 

and 

Ep = E^ = 5.83 g'^'^/L . (66) 

The lowest excited energy levels are the states of the ground state ro- 
tational band K = Q : m{2'^) = 0.35 and ?7i(4+) = 1.18, the states of the 
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p-vibrational band Up = 1 : m(0+) = 1.24, m(2+) = 1.59 and the state in the 
band Up = 2 : r7i(0+) = 2.5 (all the masses are given in units of 1/L). 

For the p vibration, we have also investigated the importance of the an- 
harmonic terms. At the lowest order in e, the collective potential is given 
by 



99 

8L^p 

or, using the appropriate rescaling p = g^/^ pjL 



V\P) = 9 P + ^TTITJ . («'') 



nfl^/'»(,^+^)i. (68) 

The cubic and quartic terms give sizeable contributions of opposite signs. If 
we approximate the collective potential by a quartic polynomial, the mass of 
the first 0"*" excited state is lowered to : m(0+) ~ 1.0/L. 

To perform a comparison with experimental data it would be necessary to 
generalize the present calculation to the SU(3) group and to include dynami- 
cal quarks. The experimental situation however appears somewhat unsettled 



5 Discussion and perspectives 

In order to compare with the approaches which consider gauge invariant 
variables constructed from the magnetic field |[l|, , it is usefull to remember 
the following relations valid in the strong coupling limit between the three 
eigenvalues &^ of the matrix BiaBja and the three A's : 

bi = 9X2 A3 , &2 = fi-Ai A3 , 63 = gXi A2 . (69) 

We have also B.y = 6 &i62&3- 

Our result for the behaviour of the wavefunction in the infrared region 
is in agreement with the conclusion of K. Johnson |I[], who found that the 
vacuum wave function should vanish at det-B = 0. The results p^ 7^ and 
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5^ 7^ at the ininiinuin cannot be obtained in a perturbative approach. 
We have found a strongly deformed minimum and the lowest energy levels 
correspond to a rotational band and a p-vibrational band. Our predictions 
for the lowest glueball states are different from those of Liischer and Miinster 
who have performed a perturbative calculation for the SU{2) gauge theory in 
a finite volume . They have choosen a basis of eigenf unctions proportional 
to exp{—^ujAiaAia) to diagonalize their effective hamiltonian. This basis 
is not adapted to describe our wave functions (|60D : a very large number of 
terms would be necessary to have a sufficient accuracy. In contrast, if we had 
obtained a minimum located at p = 0, a gaussian in the cartesian coordinates 
would have been an acceptable ansatz 0. Because the minimum is located at 
p 7^ 0, the correct ansatz is a gaussian in the "curvilinear coordinates" 6p, 6e 
and 7], which cannot be expressed as a gaussian in the cartesian coordinates 

Let us make some remarks about the derivative expansion. This is an 
expansion in powers of Xj g^l"^ . It is expected to be valid when the derivatives 
of gauge invariant quantities are small compared to some scale, for instance 
|<9iAn| ^ p^. In lowest order in the strong coupling approximation, there 
is no propagation and the wavefunctional reduces to a function of gauge 
invariant quantities constant in space. The next step is to investigate the 
coupling between the sites and its effect on the vacuum state properties and 
the glueball spectrum. 

For color singlet states of zero angular momentum, the first term involving 
derivatives in the kinetic energy is : 



^1 = - A I E ^n sL (v./") ^ A„ sL (v.f ) 

"^9 [n,m,k "^fc 



_5_ 



(70) 



where 



ok _ ^knm /-..N 

nm \2 \2 ' \ / 
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and /" are the column vectors of the matrix /. For the potential energy, 
the terms containing derivatives are exactly given by the last term in eq. 
(P). In order to check the accuracy of our variational Ansatz, it will be 
usefull to compare our results with other SU(2) calculations : finite volume 
results 0, lattice Monte Carlo calculations [§ and analytical strong coupling 
expansions . This comparison will also provide a usefull guide to construct 
variational Ansatze for negative parity states which are not included in our 
present variational space. Indeed the rescaling we have performed implies 
trial wave functions of the form \l/(Af, Xl, A3). 
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